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Dynamics near the Surface Reconstruction of W(IOO) 
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Using Brownian molecular dynamics simulation, we study the surface dynamics near the recon- 
struction transition of W(IOO) via a model Hamiltonian. Results for the softening and broadening 
of the surface phonon spectrum near the transition are compared with previous calculations and 
with He atom scattering data. From the critical behavior of the central peak in the dynamical 
structure factor, we also estimate the exponent of the power law anomaly for adatom diffusion near 
the transition temperature. 
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The (1x1) to c(2x2) structural transition on W(IOO) 
has been one of the best characterized and studied sur- 
face phase transitionadu. Besides the obvious interest in 
surface science as an important surface structural transi- 
tion, it provides a physical realization of the universality 
class of two dimensioHal XY model with four-fold sym- 
metry breaking fieldffl'Q. Moreover, the strength of the 
symmetry breaking field can be tuned by varying the 
coverage of adsorbates such as Hydrogen atoms. In par- 
ticular, the field can be made to vanish at a particular 
coverage to realize the pure two-dimensional XY modelu. 
As such, the study and understanding of this transition 
has significance far beyond this system itself. It is of 
importance in the general field of statics and dynamics 
of two dimensional phase transitions. The equilibrium 
properties of the W(IOO) transition are well understood 
through extensive experimental and theoretical investiga- 
tions. On the other hand, the critical dynamics for this 
transition has received relatively little attention. Exper- 
imentally, there exists only one study which nrobes the 
softening and damping of the phonon modesEJ near the 
critical temperature T^. However, the details of the crit- 
ical dynamics such as the emergence of the central peak 
and the dynamical critical exponents have not been inves- 
tigated. Previous theoretical studies of the dynamics for 
this transition consist of a finite slab geometry Molecular 
Dynamics (MD) simulation studyH and-an analytic con- 
tinued fraction (CF) expansion methodQ. Both of these 
approaches have limitations near Tc either because of the 
small size limit in MD or because of the truncation of the 
CF expansion after a few steps. In particular, although 



the softening and overdamping of the vibrational modes 
were confirmed in these studies, both these studies were 
not able to determine the dynamical and static critical 
exponents which characterize the critical region quanti- 
tatively. 

Recently, there are renewed interests in studying the 
critical dynamics of the W(IOO) surface reconstruction 
using different experimental probes. One approach that 
has been proposed is a time resolved laser and LEED 
spectroscopy which consists of an initial laser pulse that 
heats up the system followed by synchronized electron 
beam pulse that probes the subsequent temporal varia- 
tion of the order parameter for the c(2x2) phaseo. This 
is a difficult experiment and has riot been attempted yet. 
More recently, Xiao and AltmanQ propose to probe di- 
rectly the temperature dependence of the diffusion con- 
stant of adsorbed atoms on W(IOO) using the degradation 
of the diffraction signal from a laser-written grating on 
the adlayer. The interest in this study arises from the 
theoretical prediction that the adatom diffusion constant 
would vanish at the transition due to-jthe coupling to the 
critical fiuctuations in the substrataj. This experiment 
is under way and the results should yield better under- 
standing of critical fluctuations near Tc for this surface. 

In view of the ongoing experimental study and the in- 
adequacies of the previous theoretical investigations, we 
have undertaken a new study of the surface dynamics 
near the W(IOO) reconstruction through the direct sim- 
ulation of the Langevin dynamics of a two-dimensional 
effective Hamiltonian which has been very successful iti 
the study of the equilibrium properties of this transitioncl. 
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The approximate CFj^tudy of the dynamics is also based 
on this HamiltonianD. Unhke the previous MD stud- 
ies, the simphcity of this effective Hamiltonian aUows us 
to use much larger sizes which is important for prob- 
ing critical fluctuation effects. We also avoid the trun- 
cation approximations used in the analytic continued 
fraction studyQ. The only drawback of using a two- 
dimensional effective Hamiltonian is that the long wave- 
length phonons, which penetrates deeply into the bulk of 
the substrate, cannot be described accurately within this 
approach. However, unlike the excitations near the crit- 
ical wavevector qo = (7r/a,7r/a), these long wavelength 
phonons are scarcely affected by the critical fluctuations 
and contributes little to the strong temperature depen- 
dence of dynamical quantities near Tc- As shown below, 
with the larger size of the present study, we are able 
to collapse all our data for the dynamical structure fac- 
tor based on a general dynamical scaling analysis. This 
allows us to extract the static and dynamical critical ex- 
ponents, as well as extending the validity of the results 
arbitrarily close to T^. 

The clean W(IOO) surface is deaeribed by an effective 
Hamiltonian in dimensionless formcl 

H — y i Ui H uj 

i 

+8H4Ui^Uiy + Ci{ui ■ Uj + Ui ■ Uk) (f) 

+C2[Uix{Ujiy - Uk'y) + Uiy{Uj>x - Uk'x)]} 

where pi is the momentum vector and Ui the displacement 
vector of a W atom from the ith site, Ri = ao{l,m), with 
I and m integers, and j = {I + l,m), k = (Z,to + 1), 
j' = {I + l,m + 1) and k' ^ {I — 1, m + 1). The dimen- 
sionless parameters used in the simulation are the same 
as for model I in the previous studyO, A = — fO, B = 40, 
Ci = 3.75, C2 = Ci/2 and = -1.85. For this choice 
of parameters, the ground state displacement is uq = 1, 
and the transition temperature has been determined pre- 
viously from Monte Carlo simulationsO to be Tc = 2.11. 
The scale factors of distance Ls and temperature Tg are 
then chosen so that uq and Tc agree with the experimen- 
tally observed values, O.2A and 230X, respectively, and 
have the values Ls = 0.2 A, Ts = 109 K. This effective 
Hamiltonian has been very successful in accounting for 
the static critical phenomena such as the detailed tem- 
perature dependence and anisotropy of the intensity of 
the LEED diffraction patternQ near Tc. This Hamilto- 
nian belongs to the universality class of the XY model 
with four-fold symmetry breaking field reopjsented by 
the Hi coupling term in the HamiltonianElllj'Ll. Thus, 
the static critical exponents are nonuniversal, depend- 
ing on the value of the marginal field H^. In the large 
coupling limit, it should correspond to the Ising univer- 
sality class when C2 — 0. It can also be related, by 
renormalization-group arguments, to the Ashkin- Teller 
model which consists of two coupled Ising models and is 
known to iiave a critical line with varying static critical 
exponentstil. However, less is known about the dynami- 
cal critical behavior. It is not known for example whether 



the dynamical critical exponent z is dependent on the 
value of the four- fold symmetry breaking field H/^. The 
same question arises for the critical line of the Ashkin- 
Teller modeltj. The time scale tg for this model Hamilto- 
nian is not adjustable and is determined by the relation 
tg — {msLl/kBTsY^^ — 0.28 ps wliere mg is the W atom 
mass. In the previous CF studyQ, it was demonstrated 
that the softening and damping of the srttface phonons 
observed in the He scattering experimentti are well de- 
scribed quantitatively without any other adjustable pa- 
rameters by the results of this model up to the point 
where the soft modes become completely overdamped. 

To simulate the dynamical behavior of the system de- 
scribed by the Hamiltonian in Eq. (|^) we|-use standard 
methods of Brownian Molecular DynamicscJ. The equa- 
tions of motion for the particle coordinates are given by 

d'^Ui dui dH 

where / is a random force with each component a satis- 
fying 

< fia{t).fji3i0) 2mkBT rjo 5{t) Sij S^p 

and rjo is a microscopic damping parameter representing 
the frictional coupling of the effective two dimensional 
system to the heat bath consisting of the other subsur- 
face layers. A N x N square lattice with periodic bound- 
ary conditions has been used in the simulations. Most 
calculations were done with A'' = 32 and A^ = 46. The 
equations of motion ^ are integrated numerically with 
a time step of dt — O.Oltg in reduced units of time tg. 
The damping coefficient parameter ijq cannot be deter- 
mined from equilibrium quantities like the other param- 
eters in the effective Hamiltonian. To determine 770, we 
compare the width of the phonon modes calculated from 
the present model in a regime far from Tc with the ex- 
perimental datali. For rjo > 0.05, the noncritical phonon 
modes have widths substantially larger than the observed 
experimental value. This sets an upper bound for rjo- In 
our present investigation, we set 770 = 0.05. Near the 
transition, the critical dynamics and damping are dom- 
inated by the interparticle interactions and the precise 
value of rjo is immaterial. 

The dynamic properties of the system can be most 
conveniently studied by examining the dynamic structure 
factor Sa(3{ci,i^) defined as 

Scp{ci,uj) = j e'^' < <(q,i)"Mq) > dt (3) 

where a, (3 represent the components of the displace- 
ment. We consider in this paper mainly the longitudi- 
nal SLiq,'-^) and transverse ^'^(g, cj) components of the 
tensor S'a/3(q, w), obtained from displacements Ua along 
the direction of the vector q and perpendicular to q, re- 
spectively. For a given wave vector q, the longitudinal 
(transverse) phonon mode shows up as a peak in S'L(q, uj) 
{St{(1,uj)), while the width of the peak corresponds to 
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the damping of the mode. In Figs. 1, 2, 3 and 4, the 
behavior of the calculated dynamic structure factor as 
a function of temperature near the critical temperature 
Tc ^ 230K are shown. As can be seen from the Figures, 
as one approaches Tc from above, both the longitudi- 
nal and transverse modes near qg = (tt/a, 7r/a) start to 
soften and broaden until they are completely overdamped 
at r '--^ 1.5 Tc- As one further approaches Tc, the dynam- 
ics is completely dominated by a strong central peak. 
These results are qualitatively similar to that obtained 
in the previous study through the analytic CF expansion 
approachH. Moreover, the peak position and width of the 
soft phonon modes in the two approaches are approxi- 
mately the same as shown in Figs. 5 and 6. Thus, the 
good agreement between the theoretical results for the 
soft modes and the experimental He atom scattering data 
as demonstrated previouslyO remains true. However, as 
shown in Fig, 6, the present results show a stronger cen- 
tral peak relative to the soft phonon mode near tlip criti- 
cal point when compared with the CF calculation^. This 
is due to the fact that a truncation of the CF after a 
few steps introduces effectively a finite size limitation, 
leading to an underestimate of the critical fluctuation 
effects. This difference in the relative strength of the 
central peak could be important in such physical quan- 
tities as the anomalous temperature dependence of the 
diffusion constant as discussed below. It also shows that 
the apparent convergence of the CF approach after a few 
steps can be deceptive, and the quantitative behavior of 
the critical dynamics needs a long expansion of the CF 
which may render the analytical approach impractical. 

We now examine the critical properties more closely 
by focusing on the central peak. It follows from standard 
scaling arguments that near Tc, the dynamical structure 
factor should satisfy the scaling form 

iV'^5(g,c^)=r+^/^5±(9^,^r) , (4) 

where q — |q — Qol, 5±(^) is a scaling function, ^ cx 
\T/Tc — is the divergent correlation length and d is 
the system dimension. We now focus on S{uj) = S{q — 
OjUj). The results for S{uj) from the present calculation 
for various uj and T near Tc are plotted in Fig. 7. Note 
that at q = qo or g = 0, the longitudinal and transverse 
modes are identical and S(uj) is isotropic. According to 
the scaling form in Eq. (||), the data in a scaling plot 
S{uj) \T/Tc - 1]" X u) \T/Tc - 1|-^ " should col- 
lapse on to a single curve above or below Tc respectively 
if the exponents u z -\- ^ and zv and the critical temper- 
ature Tc are chosen correctly. This scaling plot is shown 
in Fig. 8 for a system size N = 32. By adjusting the 
parameters for the exponents and Tc, we have verified 
that the data for S{uj) indeed obeys the scahng hypoth- 
esis. We obtain the values Tc = 2.07(4), vz = 1.9(4) and 
vz -\- ^ = 3.8(4). The same analysis was performed for a 
larger system size iV = 46 with the same results indicat- 
ing that finite-size corrections are negligible in the range 
of frequency and temperatures studied. The resulting 
value of 7 = 1.9(4) is consistent with previous estimate 
from equilibrium simulations which foundlj 7 = 1-5(2) 



and V = 0.83(1). Combining this value of ly with the 
present results from dynamic scaling, we get an estimate 
for the dynamic critical exponent z = 2.3(4). It is in- 
teresting to note that according to the dynamical scaling 
arguments, S{uj = 0) diverges as |T/Tc - Ij-^" ^+'') as 
one approaches Tc. In the analytic CF approach how- 
ever, a truncation of the continued fraction for S{lu = 0) 
at the second stage yields an approximate divergent be- 
havior of |T/Tc - 1|-2t for S{uj = 0). For the present 
model, these two exponents ly z + j and 27 are indistin- 
guishable within our numerical error. Another interest- 
ing observation is that the value of the dynamical critical 
exponent z determined here is close to the correspond- 
ing value for the Ising model which lies in the range 2.08 
to 2.24 as obtained, by numerical and renormalization- 
group calculationstJ. It raises the interesting question 
of whether the dynamic exponent can be universal even 
when the static critical exponents are known to be de- 
pendent on the value of the marginal field H4. Further 
work on dynamics of the XY model with four-fold sym- 
metry breaking using more efficient Mopte-Carlo meth- 
ods for larger size systems are underwayfij. They should 
help to clarify more the issue of dynamic universality 
class for this model. However, for thp-j-elated Ashkin- 
Teller model, recent numerical resultgl3 show that the 
dynamical exponent z does not vary significantly along 
the nonuniversal critical line compared with the error es- 
timates and is also close to the results for the Ising model. 
Thus, we expect that the present estimate of z may be 
a good approximation for our model with a weak depen- 
dence on the parameters of the model. 

While the soft modes and the central peak can be di- 
rectly probed experimentally with He scattering studiesEJ, 
it also influences indirectly physical quantities such as 
adatom diffusion constant which depends strongly on the 
coupling of the adatom to the substrate excitations. As 
shown previouslyB, the friction rj acting on the adatom 
due to the coupling to the substrate phonon excitations 
is related to the zero-frequency limit of the dynamical 
structure factor, 

where W{q) denotes the Fourier coefficients of the cou- 
pling of the adatom to, the substrate. Thus, adatom 
diffusion measurementsQ will probe the dynamic corre- 
lation function S'(g, w = 0) which diverges as the tem- 
perature approaches Tc because of the central peak exci- 
tations. According to the general scaling relation given in 
Eq. (^, the friction coefficient rj should diverge near the 
transition as |T/Tc — 1]~^, with the critical exponentt£l 
X — v{z — 2) -|- 7. The difference of the exponent for 77 
from that of 5'(a; = 0) arises from the integration over the 
wave vectors for the friction coefficient. Since D cx l/ry, 
this leads to the prediction of a vanishing diffusion con- 
stant as one approaches Tc, with a power-law anomaly 
D cx |T/Tc — ll"^. The divergent behavior of the fric- 
tion parameter near Tc is not just specific to this system. 
Diffusion anomalies near a surface preroughening transi- 
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tion has also been predictedtZl based pa similar considera- 
tions. Finally, we have shown ear lierta that, the effective 
pinning strength of flux lattice in disordered type 11 su- 
perconductors should diverge at the flux-latticc melting 
transition due to the coupling of the pinning center to the 
critical fluctuations. This could be the possible origin for 
the "peak effect" observed in high-Tc supercpiiductors 
in a magnetic field. Recent experimental dataEfl show a 
precise correlation of the "peak effect" with the vortex 
lattice melting transition, lending support to this picture. 

In conclusion, we have performed a Langevin simu- 
lation study of the dynamics of the W(IOO) surface re- 
construction near the transition point. We are able to 
avoid the severe size limitations of the previous molec- 
ular dynamical studies by studying an effective two di- 
mensional Hamiltonian and analyse the results through 
dynamical scaling theory. We have determined both the 
static and dynamical critical exponents for this transi- 
tion which belongs to the universality class of XY model 
with four fold anisotropy. Finally, we have related the 
central peak of the phonon spectrum to the friction ex- 
perienced by an adatom on this surface. The adatom dif- 
fusion constant on this surface is predicted to vanish at 
the transition point due to the divergent frictional damp- 
ing. The comparison of our new results with forthcoming 
adatom diffusion data would provide a detailed picture 
and further understanding of the critical dynamics near 
this intriguing surface reconstruction. Furthermore, the 
divergent nature of the friction near the transition is not 
just specific to this system but should manifest itself near 
other continuous phase transitions such as the surface 
preroughening transition and the vortex lattice melting 
transition in type II superconductors. 
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98/00977-2 (M.R.B). 
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FIG. 1. Longitudinal structure factor ^(q, w) for 
q = 0.6qc, at different temperatures. 



FIG. 2. Transverse structure factor St{ci,ijj) for q = 0.6qo 
at different temperatures. 
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FIG. 3. Longitudinal structure factor ^(q, w) for FIG. 4. Transverse structure factor St (q, w) for q = 0.8qo 

q = O.Sqc, at different temperatures. at different temperatures. 
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FIG. 5. Longitudinal structure factor Sl{ci,i^) far from Tc FIG. 6. Longitudinal structure factor Sl{ci,ui) close to Tc 

for q — O.Sqo ff obtained from the present study and the CF for q — 0.6qo ff obtained from the present study and the CF 
approximatiorjj. approximationu. 
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FIG. 7. Structure factor S{u!) at q = qo as a function of 
dimensionless temperature T and frequency ui near Tc . The 
units for T and u> are Tg and tj^ as described in the text. 



FIG. 8. Scaling plot of SMIT/Tc-ir^+^xwIT/Tc-ir"^ 
from the data in Fig. 7 obtained by adjusting the parame- 
ters Tc, vz and vz + 7. The upper set of collapsing data 
corresponds to T > Tc and the lower one to T < Tc. 
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